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Recent numerical work by Wilson, Mathews, and Mar- 
ronetti [J. R. Wilson, G. J. Mathews and P. Marronetti, Phys. 
Rev. D 54, 1317 (1996)] on the coalescence of massive binary 
neutron stars shows a striking instability as the stars come 
close together: Each star's central density increases by an 
amount proportional to 1/ (orbital radius). This overwhelms 
any stabilizing effects of tidal coupling [which are proportional 
to 1/ (orbital radius) 6 ] and causes the stars to collapse be- 
fore they merge. Since the claimed increase of density scales 
with the stars' mass, it should also show up in a perturba- 
tion limit where a point particle of mass [s, orbits a neutron 
star. We prove analytically that this does not happen; the 
neutron star's central density is unaffected by the compan- 
ion's presence to linear order in fi/R. We show, further, that 
the density increase observed by Wilson et. al. could arise as 
a consequence of not faithfully maintaining boundary condi- 
tions. 

PACS numbers: 97.80.Fk, 04.25.Dm, 04.40.Dg, 97.60.Jd 



Wilson, Mathews, and Marronetti (WMM) jy have 
proposed a method of approximating the fully General 
Relativistic analysis of binary neutron star coalescence. 
The essence of their scheme is to choose a simple form of 
the spacetimc metric (one in which the spatial three slices 
are conformally flat), and solve the constraint equations 
of General Relativity (GR) for some initial matter con- 
figuration. They evolve only the fluid equations forward 
in time until the fluid reaches a quasi-equilibrium con- 
figuration, then solve the constraint equations again for 
the new matter configuration and iterate until a quasi- 
equilibrium solution to the combined Einstein-fluid equa- 
tions is found. Their method makes 3-dimensional sim- 
ulations of such systems more tractable by reducing the 
computational requirements. 

These simulations yield an extremely surprising result: 
neutron stars that are close to the maximum allowed 
mass are "crushed" into black holes long before the neu- 
tron stars coalesce. They claim that the origin of this 
effect is a non-linear gravitational interaction due to the 
companion's presence that strengthens the gravitational 
potential of each star. Consider a binary star system — 
star- A has mass Ma, and star-B has mass Mb- WMM 
claim that non-linear interactions cause the potential at 
star- A to be increased by a term that scales as Mb/R 
(where R is the orbital separation). This, in turn, in- 
creases the internal energy and density of star-A by terms 
that scale as Mb /R. If star-A happens to be marginally 
stable in isolation, the effect is sufficient to push it over 



the edge, causing a catastrophic collapse to a black hole. 

The scaling law claimed by WMM is precisely what one 
would expect if this effect were due to a post-l-Newtonian 
enhancement of the gravitational interaction. Motivated 
by this observation, Wiseman has recently done a careful 
analysis of the effect that a binary companion has on a 
fluid star, using the first post-Newtonian approximation 
to GR H ; he finds no change to either the central energy 
density or the angle averaged proper radius of the star 
at this order. Wiseman's calculation does not rule out 
completely a star-crushing effect, but does show that it 
is not evident at post-l-Newtonian order in GR. 

Suppose for a moment that the WMM effect is a prop- 
erty of neutron star binaries in GR, and that it scales as 
Mb/R at star-A. Clearly, it should also be apparent in 
the limit that we shrink star-B down to a point particle 
of mass /i <C Ma = M. In this limit, the exact solution 
of the Einstein field equations describing a binary neu- 
tron star system can be approximated by a perturbative 
expansion in p,/ R about the solution for an isolated star. 
We write the metric as g a p — g^g + eh a p + 0(e 2 ), where 
the superscript indicates the background metric, and 
we have introduced an order counting parameter e with 
the formal value unity. Quantities multiplied by e scale 
linearly with p/ R, quantities multiplied by e 2 scale with 
(p/R) 2 , etc. In what follows, we ruthlessly discard all 
terms of order e 2 , constructing an argument that is valid 
only to linear order in fj,/R. 

The neutron star material is considered to be perfect 
fluid with stress-tensor 



T a /3 = Pg a p + (P + p)u a Uf} . 



(1) 



This must be supplemented with an equation of state 
relating the pressure P and the energy density p. For 
concreteness, we assume a poly tropic form P = Kn r 
where K and T are constants, and n is the fluid's baryon 
density. The energy density p is directly related to n 
by the first law of thermodynamics; see Eqs. (3.2.6-7) of 
Ref. 0. 

We take the background spacetime to be that of an 
isolated, spherical star with the line element 



ds 2 = -e 2 *^dt 2 



dr 2 



[1 - 2m(r)/r] 



*dn 2 



(2) 



Here m(r) is the gravitational mass inside a sphere of 



radius r, and dfl 2 



de 2 



sin 2 9 dip 2 . The combined 



Einstein-perfect fluid equations, generally referred to as 



1 



the Oppenheimer-Volkoff (OV) equations (see for exam- 
ple Chapter 23 of Ref. j|) are solved by demanding regu- 
larity of the origin [m(0) = 0] and fixing the value of the 
central baryon density n c . The radius of the star R'g is 
the coordinate radius at which the baryon density n° be- 
comes zero. The central density n c uniquely determines 
the total mass M = m(-R|§), and baryon mass My, 
this statement is equivalent to Theorem 7 of Ref. || . 

One can show using the polytropic equation of state 
and the OV equation for the pressure, 



dP° 
dr 



(p° + P°)[m(r) + 47rr 3 P°] 



r\r 



2m(r)] 



that 



dn° 
dr 



(3) 



(4) 



Finally, the background geometry outside the star is de- 
scribed by the Schwarzschild solution with m = M and 
exp(2$) = 1 - 2M/r in Eq. (|). 

The perturbing source is a single point particle of 
proper mass /x in a circular orbit at radius R. It is de- 
scribed by the stress-energy tensor 



T 



a0 _ e^f_ 
' R 2 i 



-S(r — R)S(cos 



fit), 



(5) 



where v a = (l-SM/R)- 1 / 2 ^, 0, 0, ft) and ft = y/M/R 3 . 
The presence of this point "star" will alter the geometry 
and disturb the material in the central star, modifying 
the description of the spacetime and matter by terms of 
order e. Linearizing G a p = %nT a p and T"' 3 ;j g = in e, 
we find that the first-order perturbation equations sep- 
arate by expanding the angular dependence in spherical 
harmonics, and the time dependence in Fourier modes. 
This is enough to address the issue of how the central 
density scales. 

Consider the expansion of the baryon mass density. It 
may be written 



n(r, 9, 0, t) = n°(r) + eV 8n lmuj (r)Yi m (9, 4>)e 



l.m.UJ 



(6) 



An immediate consequence of Eq. (^) is that 5ni muJ (0) = 
for I > 1: if it were non-zero, the density would be 
multi- valued at r = 0. Thus, only the monopole could 
affect the central density if the center of the perturbed 
star were to remain at the origin. In reality, the star's 
center star will not be at the coordinate origin: the or- 
biting body will move it to some point in the orbital 
plane. However, the magnitude of this shift must be of 
the same order as the perturbation itself: r ccnt = e£(t) 
for some function £(t). Now evaluate the density at the 
star's center with a Taylor expansion: 



n dn° 
n CO nt= n (0) + e£ — — 
dr 



r=0 



Sni moj (0) + e£ 



dSnir, 



dr 



r=0. 



Ylr. 



(7) 



As we have already shown, Sni mLJ (0) = 0, except possibly 
for I = 0, while dn°/dr -> as r -> by Eq. (g). Thus, 
the baryon density at the center of mass is given by 



n°(0) + efe 00 o(0) + 0(e 2 ) 



(8) 



Only the monopole can produce changes in the central 
density which scale linearly in \xjR. 

It is straightforward to solve for the I = m = correc- 
tions to the mctricoutsidc the fluid. Define the function 





H(r) = { 2 f i (1-2M/R) 
r (1-3M/R) 1 / 2 



r <R 
r>R 



(9) 



then htt = H(r), h rr = H(r)/(1 - 2M/r) 2 . We have 
set hee — = using the first order gauge freedom 
available for the monopole. As the point particle spirals 
into the star, some of its total mass-energy (its contri- 
bution to the total mass as measured at infinity) is radi- 
ated away, though its locally measured mass (rest mass) 
is conserved. The multiplicative factor in the above ex- 
pression correctly accounts for that radiation loss. Notice 
that there is no monopole contribution to the metric in- 
side the orbital radius of the particle — Keplerian orbits 
inside the orbital radius measure only the mass of the 
unperturbed neutron star at monopole order. 

Is it possible for the monopole part of the perturbation 
to rearrange the fluid in the star, but leave its total grav- 
itational mass unchanged? The answer is unequivocally 
no. A monopole perturbation is spherically symmetric 
and can only take one spherical solution into another. 
However, when the equation of state is fixed, all spherical 
solutions are parameterized by the gravitational mass — 
for each value of the gravitational mass M there exists 
a unique spherical configuration of the star. Spherical 
solutions therefore exhibit a one-to-one correspondence 
between the gravitational mass and the central density of 
the star. Since the gravitational mass M is unchanged at 
monopole order, the central density cannot be affected ei- 
ther. Hence, the central energy density of a neutron star 
is unaffected by a binary companion at order fj,/ R. There 
is no crushing effect which scales linearly with \xjR. 

Interestingly, it is easily shown that incorrectly im- 
posing boundary conditions can lead to an increase in 
central density at order fi/R. If the total gravitational 
mass of the star and particle is held fixed in a sequence 
of quasi-equilibrium solutions (ignoring the gravitational 
radiation that causes the orbital radius to shrink), and 
the particle's locally measured mass (rest mass) is held 
fixed, then obviously the star's total mass and baryon 
mass must go up by an amount of order \iM / R, contrary 
to how a real binary would behave. This mass increase 
will drive the central density up by a fractional amount 
of order /i/R, which is what the WMM simulations show. 
We have no evidence that this is what actually happens 
in the WMM simulations; it merely illustrates one way 
in which the observed density increase could arise. 
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